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Abstract 

We  studied  the  contribution  of  internal  energy  and  entropy  change  to  cathode  voltage  assuming  that  the  cathode  material  Li,-Mn204  was 
completely  ionic  and  only  the  Coulomb  potential  was  effective.  We  calculated  the  cathode  voltage  of  an  ideal  spinel  structure  LitMn204 
(Fd3m)  (Mn-O:  2  A,  Li-O:  1.732  A,  O-O:  2.828  A).  From  this  calculation,  we  found  that  LivMn204  has  a  multiphase  co-existence  for 
0  <  x  <  1  with  a  very  high  voltage  (8.7  V  versus  Li),  and  has  a  negative  voltage  when  1  <  x.  If  we  suppose  that  inserted  Li+  and  electron 
have  40%  of  the  formal  charge,  LivMn204  has  a  voltage  from  4.55  to  4.79  V  for  0  <  x  <  1,  and  has  a  voltage  of  1.75  V  when  1  <  x. 
©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 


2.  Basic  theory  and  model  [6,7] 


Recently,  there  have  been  several  studies  predicting  the 
cell  voltage  of  cathode  active  materials  in  lithium  ion  cells 
based  on  the  first  principles  quantum  chemistry  [  1— 5] . 
However  at  present,  it  is  not  easy  to  calculate  non-averaged 
voltage  using  the  first  principles  quantum  chemistry  because 
of  long  calculation  time. 

We  have  been  studying  the  possibility  of  predicting  the 
voltage  from  the  Coulomb  potential.  The  cell  voltage  is 
described  thermodynamically.  We  considered  the  contribu¬ 
tion  of  the  internal  energy  change  to  the  cell  voltage  (E 
versus  Li/Li+)  assuming  that  the  cathode  material  was 
completely  ionic  and  only  the  Coulomb  potential  was 
effective  in  terms  of  changing  the  internal  energy.  We 
calculated  the  cathode  voltage  of  an  ideal  structure  Li,M02 
( R3m )  (Li-O:  2  A,  Ni-O:  2  A)  [6,7],  which  is  a  model 
compound  of  LiTNi02.  From  this  calculation,  we  found  that 
Li  VM02  has  a  two-phase  co-existence  for  0  <  x  <  0.25,  and 
one  phase  for  0.25  <  x  <  1. 

In  this  paper,  we  calculate  the  cathode  voltage  of  an  ideal 
structure  spinel  LiYMn204  ( Fd3m )  using  the  same  method 
employed  for  the  calculation  [6,7]  of  LivM02  (R3m). 
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We  assume  that  the  cathode  active  material  is  an  ionic 
crystal  and  that  the  charge  is  at  the  center  of  each  atom.  We 
ignored  any  forces  between  atoms  other  than  the  Coulomb 
force. 

Let  >-jj  (A)  be  the  distance  from  a  certain  atom  k  to  the 
surrounding  atoms  of  species  i  in  an  ionic  crystal,  and  Zi  be 
the  charge  number  of  species  i.  Then  the  Coulomb  potential 
at  atom  k  is  expressed  by  the  following  Eq.  (1). 


^  =  ^E-=14-3"E- 

47t£04fr;i/  y-fru 


(1) 


where  (Pk  (eV)  is  the  Coulomb  potential  at  atom  k,  e  the 
elementary  charge,  and  s0  is  the  permittivity  of  free  space. 
From  thermodynamics,  we  have 


(2) 


where  E  is  cell  voltage  (V),  AGr  the  reaction  Gibbs  function 
(J/mol),  n  the  electrons  per  molecule  oxidized  or  reduced, 
and  F  is  the  Faraday  constant  (C/equiv).  From  Eqs.  (1)  and 
(2),  we  have 
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where  C  is  the  reaction  coordinate  (mol),  P  the  pressure  (Pa), 
T  the  temperature  (K),  U  the  internal  energy  (J),  V  the 
volume  (m3),  and  S  is  the  entropy  (J/K).  Let  t/a  (eV,  energy 
change  for  one  atom  Li  insertion)  be  the  internal  energy 
change  caused  by  the  insertion  of  one  Li.  Then  we  have,  for 
the  first  term,  Eml  =  —UJne.  We  will  ignore  the  second 
term,  because  the  value  is  very  small.  The  third  term  (£ent) 
is  the  entropy  term:  Z?ent  =  (RT/nF)  In  ((1  —y)/y),  where 
y  is  the  probability  of  Li  occupation  in  cathode  Li  sites.  The 
electron  arrangement  does  not  contribute  to  the  entropy 
term  because  we  assumed  that  electron  is  inserted  into 
the  nearest  Mn  from  the  inserted  Li+. 

We  considered  the  energy  cycle  for  the  reaction 

(5£Li(crystal)  — >  (5£Li+(g)  + 
and 


(5£Li+(g)  +  <5£e  +  Li,Mn204  -►  Liv+^Mn204 


Let  Ad  be  the  ionization  potential  (eV)  of  Li  metal,  and  (I>t  l 
be  the  Coulomb  potential  at  a  vacant  Li+  insertion  site.  The 
increase  in  internal  energy  caused  by  the  insertion  of  an 
electron  is  assumed  to  be  —  <Pmr  +  /mh,  where  <PMn  is  the 
Coulomb  potential  of  the  electron  insertion  site  and  /Mn  is 
the  LUMO  energy  of  an  isolated  Mn  atom  in  a  vacuum. 
When  <PMn  is  calculated  for  a  cathode  before  the  Li+ 
insertion,  we  have  to  add  the  contribution  of  the  inserted 
Li+.  Thus,  we  have  to  add  1 4.399//^r  i  _Mn  to  <PM,  where 
^u-Mn  is  the  distance  between  the  inserted  Li+  and  the 
electron.  Finally  we  have  the  following  equation  where  the 
sign  after  /Li  is  +  if  there  is  already  Li  at  the  center  of 
k  (extraction  of  Li),  and  —  if  there  is  no  Li  at  the  center  of  k 
(insertion  of  Li).  With  Li  ion  cells,  n  =  1. 
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3.  Calculation  of  cathode  voltage  of  LixMn204 

In  this  paper,  we  assume  an  ideal  spinel  LiMn204  which  has 
the  same  nearest  length  of  2 r  (A)  for  0-0  or  Mn-Mn.  In  this 
ideal  material,  oxygen  atoms  form  regular  tetrahedrons  linked 
at  their  edges  to  form  oxygen  layers,  and  MnOe  and  Li04  form 
regular  octahedral  and  tetrahedral  structures,  respectively.  We 
used  r  =  2°  5 .  Then,  the  oxygen  layer  distance  of  the  ideal 
structure  is  2.309  (A),  and  that  of  the  experimental  value  is 
2.38  (A)  [8],  The  distance  from  the  experiment  is  1.03  times 
greater  than  that  of  the  ideal  structure. 

The  Coulomb  potential  at  atom  k  (Eq.  (1))  is  calculated 
using  VIP  BASIC  (Mainstay)  for  Power  Macintosh. 

We  obtained  an  /Li  value  of  —7.107  eV  (685.78  kJ/mol) 
from  the  enthalpy  change  for  Li(cr)  — >  Li+(g)  +  e. 

Using  the  DV-Xa  molecular  orbital  method  (SCAT)  [9], 
we  determined  that  the  HOMO  energy  of  [MnOg]-1  cluster 
was  —0.40627  eV,  where  the  nominal  valences  of  Mn  and  O 
are  3  and  —2/3,  respectively.  The  contribution  of  the  Cou¬ 
lomb  potential  of  oxygen  is  —28.798  eV.  Therefore, 
/Mn  =  -0.40627  -  28.798  =  -29.204  eV. 

The  calculated  voltages  (Eint)  are  summarized  in  Table  1. 
When  all  the  Li+  atoms  are  moved  from  tetrahedral  sites  to 
octahedral  sites,  the  spinel  structure  changes  to  a  rock  salt 
structure. 

After  the  addition  of  entropy  term  in  Eq.  (4),  the  voltage 
profile  is  calculated  using  “Mathematica  (Wolfram)”,  and 
shown  in  Fig.  1.  The  Li  sites,  where  Li  can  be  inserted  or 
extracted,  are  C4-2  at  Z  =  3  + 1  /4,  B4-2  at  Z  =  1  + 1  /4,  A4- 1  at 
Z  =  5  +  1  /4  for  0  <  x  <  0.25,  A4-2  at  Z  =  3  +  l/4,C4-lat 
Z  =  1  +  1  / 4,  B4-2  at Z  =  5  +  l/4for0.25  <  x  <  0.5.C4-1 
at  Z  =  3  +  1/4,  B4-1  at  Z  =  1  +  1/4,  A4-2  at  Z  =  5  +  1/4 
for  0.5  <  x  <  0.75,  A4-1  at  Z  =  3  +  1/4,  C4-2  at  Z  =  1  + 
1/4,  B4-1  at  Z  =  5  +  1/4  for  0.75  <  x  <  1  (Fig.  2).  In 
each  area,  Li+  ions  are  assumed  to  distribute  at  random.  This 
means  that  all  the  sites  are  occupied  by  Li v+  and  e  v_,  where  y 
is  the  site  occupancy  in  each  area.  Then,  for  Eml,  we  have  a 
straight  line  connecting  two  voltages  (Zqnt)  at  each  end  point 


Table  1 

The  calculated  voltage  of  Li^Mi^CVLi  cell 


x  in  LixMn204 

Li  site 

Mn  charge 

Ein,  (V) 

Voltage  from 
Sint  (V) 

Ebmodint(V) 

Voltage  from 

^bmodint  (V) 

Li-Mn  layer 

Mn  charge 

Spinel 

0 

Insertion 

Tetrahedral 

+4 

+4 

6.93 

8.74 

3.81 

4.79 

0.25 

Extraction 

Tetrahedral 

+3,  +4 

+4 

9.27 

5.57 

Insertion 

9.09 

4.15 

0.5 

Extraction 

Tetrahedral 

+3 

+4 

9.65 

5.63 

Insertion 

7.00 

7.91 

3.82 

4.66 

0.75 

Extraction 

Tetrahedral 

+3 

+3,  +4 

8.81 

5.50 

Insertion 

6.73 

7.19 

3.78 

4.55 

1 

Extraction 

Tetrahedral 

+3.5 

+3.5 

7.64 

5.31 

Rock  salt 

0 

Insertion 

Octahedral 

+4 

+4 

2.03 

5.18 

- 

1.75 

1 

Insertion 

Octahedral 

+3.5 

+3.5 

5.18 

5.18 

- 

2 

Extraction 

Octahedral 

+3 

+3 

4.42 

4.42 

- 
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x  in  LixMn2C>4 

Fig.  1 .  The  calculated  voltage  of  LiA.Mn204. 


of  the  area.  We  obtained  the  following  equation  for  spinel 

E  =  £f+X,y-^to(T^;)  (5) 

where 


£?  =  6.93  (V), 

Ki  =  2.34  (V) 

and 

-V  =  025  ’ 

for 

£0  =  9.09  (V), 

K2  =  0.56  (V) 

and 

(x  -  0.25) 

for 

}  “  0.25  ’ 

£0  =  7.00  (V), 

Ki  =  1.80  (V) 

and 

(x  —  0.5) 

}  “  0.25  ’ 

for 

£0  =  6.73  (V), 

K4  =  0.91  (V) 

and 

(x  —  0.75) 
y  —  n  ’ 

for 

From  Eqs.  (2)  and  (5),  the  corresponding  Gibbs  free  energy  is 
G  =  -  ^  |  £?; y  +  ^-y2  -  ^  [y  In  y  +  ( 1  -y)  ln(  1  -y) + Q]  j 
where 


The  relation  is  shown  in  Fig.  3.  There  are  three  tangential 
lines  (dotted  lines  in  Fig.  1),  which  show  phase  separation  at 
0  <  x  <  0.5,  0.5  <  x  <  0.75,  and  0.75  <  x  <  1.  From 
Eq.  (2),  the  potential  of  the  phase  separation  area  is  —(slope 


0  <  x  <  0.25 
0.25  <  x  <  0.50 
0.50  <  x  <  0.75 
0.75  <  x  <  1, 


of  the  tangential  line).  Therefore,  the  voltage  of  bold 
curve  shown  in  Fig.  1  is  the  final  result  from  this  calcula¬ 
tion.  Unfortunately,  the  state  at  x  =  1  is  the  most  stable 
phase  and  we  could  not  observe  the  discharge  curve  for 
1  <  x  <  2. 


Ci  =  0  and  y  = 


0.25’ 
and  y  = 


(x  -  0.25) 


Ki 
2 

C4  =  C3  -  (  E'i  +  ^ 


and  y  = 
and  y  = 


0.25  ’ 

(x  —  0.5) 


0.25  ’ 

(x  -  0.75) 
025 


for  0  <  x  <  0.25 
for  0.25  <  x  <  0.50 
for  0.50  <  x  <  0.75 
for  0.75  <  x  <  1 


If  Li+  ordering  does  not  occur  for  rock  salt,  G  is 


G  =  —F<  2.03x 


2.39x2  fRT\ 


—  (  —  J  [x  ln(x)  +  (2  —  x)  ln(2  —  x)  —  2  In  2] 


0  0.5  1.0  1.5  2.0 

x  in  LixMn204 


0  0.5  1.0  1.5  2.0 

x  in  LixMn2C>4 


Fig.  3.  The  Gibbs  free  energy  of  Li  tMn20 
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Our  calculated  voltage  is  too  high  compared  with  the 
experimental  value  (about  4.5  V  (0  <  x  <  1 )  and  3  V 
(1  <  x  <  2))  [8].  We  used  40%  of  the  formal  charge  of 
inserted  Li+  and  electron  for  spinel  to  improve  the  calcu¬ 
lated  value.  The  improved  voltage  (£Wodint)  is  shown  in 
Table  1.  We  used  the  equation;  £bmodint  =  0.16£)nt  +  2.07 
(for  insertion)  and  -Ebmodint  =  0.1 6E{nt  +  4.09  (for  extrac¬ 
tion).  We  assumed  100%  of  the  formal  charge  of  inserted 
Li+  and  electron  for  rock  salt  because  good  fitting  of  LiNiC>2 
voltage  using  this  assumption.  We  obtained  G/F(spinel)  = 
— 4.70eVatx  =  1,  and  G/F(rock salt)  =  — 6.45eVatx=  2. 
Therefore,  there  is  a  two-phase  region  at  1  <  x  <  2,  and  the 
voltage  is  1.75  (V).  The  result  is  also  shown  in  Table  1. 
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